Unitary Fermi Gas: Scaling Symmetries and Exact Map 
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The weakly bound system of dimers formed in a cold Fermi gas at infinite scattering length 
reveals an underlying SU{1, 1) symmetry, which facilitates an exact map from the interacting to the 
non-interacting system. A systematic treatment of this symmetry for the ground state, is shown 
to result in a shift of energy proportional to the scaling exponent. For the excited states, novel 
breathing modes at integral values of the harmonic frequency uj are predicted in one-dimension and 
consequences for the realization of exclusion statistics are pointed out. 

PACS numbers: 



Universality of the dilute Fermi gas at unitarity has 
garnered considerable attention in the recent litera- 
ture |l| . The unitarity regime characterized by a diverg- 
ing scattering length a, is the cross section reaching a 
maximum value 47r/fc^, where k is the relative momen- 
tum of the scattering particles. Using Feshbach reso- 
nance, the infinite scattering length is achieved when the 
energy of a pair of scattering atoms is tuned close to that 
of a molecular bound state by an external magnetic field, 
leading to a substantial resonance scattering. A hallmark 
feature of this system is that the pairing correlations have 
a range shorter than the de Broglie wavelength, leading 
to significantly novel behavior 0, Q . The unitary regime 
is described by weakly interacting particles, which pro- 
vides an ideal ground to study physics at the interface 
of weak coupling BCS superfluidity and BEC. However, 
the lack of a small parameter poses a challenging problem 
in understanding the structure of many-body wave func- 
tions and the associated bound states. The only available 
length scale in the ground state is inter particle spacing, 
ri~i/3^ where n is the density, since, a must drop out 
from all physical observables. Thus, universality requires 
that for both bosons and fermions, the only relevant en- 
ergy scale is the Fermi energy Ef. The energy Eu of 
the unitary fermions has to then scale as £,Ef, where the 
parameter ^ does not depend on the specific form of the 
short range potential and is same for all potentials which 
have a zero energy S-wave bound state. Determining this 
scaling parameter ^ is central to an understanding of the 
degenerate Fermi gas and also its superfluid transition 
temperature. From the theoretical view point, a careful 
inspection [il] of the short distance behavior of the two- 
body wave function ip, in the unitary regime, shows a 
scaling behavior Q, 
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a shift in the ground state energy by a factor 

= (7 + dN/2)hw . (2) 

To study this scaling regime, it is instructive to con- 
sider the problem of short-range interacting particles in 
d-dimensions and the general Hamiltonian 
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H = y ^ 



(3) 



which satisfies the condition ([T]) and the Schrodinger 

equation -^^^^^ V^* + E,<j ^(^1 " r,)^ = 
without an external potential. Note that the two-body 
potential V depends only on the inter particle distance 
vector Ti—Tj and is characterized by the s-wave scattering 
length a. The precise form of this potential is unimpor- 
tant in the unitary regime. The scaling properties of this 
system are as follows 
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l'(r) ^ A''/2^(Ar) 
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where 7 is the scaling exponent. Amusingly, this leads to 



where we concentrate on potentials which have the scal- 
ing law y(Ar) = V{v)l}?. As pointed out in Ref. [13], 
in additional to a type potential in d-dimensions, 
it is also possible to have V{r — r') — (5^(r — r') in two 
dimensions. It should be mentioned that a distinctive 
role of c? = 4 has emerged in recent literature due to the 
1 /r^ behavior of the two-body wave function, leading to 
a logarithmic divergence in the normalization J^. Both 
numerical Q and e-expansion around d = 4 indicate, that 
the unitary regime can be well described by weakly in- 
teracting fermions and bosons t7|. Here, the potential V 
possesses a degree —2 characteristic of l/r^, as the man- 
ifestation of scale invariance can only allow a 1 /r^ type 
interaction between pairs. 

In a harmonic trap, the weakly interacting particles 
can show several new features. The interactions of the 
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type given in ([3]), can capture the dimerization effect 
accurately. The formation of molecules, which subse- 
quently condense to a BEC, necessitates the search of a 
zero-energy bound state, leading to a dimerization of the 
fermionic gas. If we consider the trapping Hamiltonian 
to be 
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then there can be additional breathing modes with uni- 
versal frequency 2a; [1, Q. Remarkably, there exists a 
hidden 50(2,1) Lorentz symmetry ^0|], the continuous 
representations of which, can explain exactly the fre- 
quency of breathing modes. The starting point for such 
an 5*0(2, 1) algebra is a suitable combination of the gen- 
erators given in Eq. (P), and (P [10]. 

In this letter, we point out that in the unitary regime 
and for a particular class of wave-functions, this underly- 
ing 50(2, 1) symmetry goes much deeper in establishing 
an exact map from the interacting to a non-interacting 
system and also naturally explains the shift in the ground 
state energy Q. The 50(2, 1) algebra in the scaling 
regime can have interesting consequences for the exis- 
tence of fractional exclusion statistics ll|. There have 



been several works in recent literature which suggest a 
number of universal features connecting interacting and 
non-interacting regimes 

[ilSIil- We wish to see which 
of these universal features are captured by the 50(2, 1) 
symmetry algebra. 

To illustrate our proposal, let us start with a general 
Hamiltonian of the type H+Htrap, having a ground state 
wave function 'i'o such that V{r) = €o + ^- Eili ^f^o, 
where, eo is the ground-state energy 1^, 3]. Let us now 
focus on a certain class of wave functions of the product 
type "^Iq = ipG, where ip captures the pairing correla- 
tions associated with the scaling properties of the dimer 
potential given in Eq. ([4]) and G governs the dynamics 
of the trap. This form of the wave function is quite suit- 
able in the dimerization regime, where we are interested 
in extracting out the features universal to the interacting 
and the non- interacting systems. It is worth mentioning 
that Ei ^ completely decouples from all the degrees of 



freedom [10|. Without loss of generality, we consider 



4' = Y[f{r 
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to satisfy the scaling behavior advocated in Eq. ([T]), 
where /(fij) denotes a homogeneous polynomial of defi- 
nite degree. For instance, if ip is of the form Y[i<j ~ 
rj\^, where i, j — 1, ■ ■ ■ , N , then the degree of the poly- 
nomial is (3N{N - l)/2. 

In order to identify the 50(2,1) or SU{1,1) genera- 
tors, it is useful to make a similarity transformation with 
respect to -0 and obtain 



(7) 



where A = ^EiliV? + Eti V,(ln ^A).V,. Let us 
note that there are two different ways of forming a 
5J7(1, 1) algebra. An inter particle potential of the 
type V{ri - fj) = ^ Ei ji'^i - ^iMVi - ^jl > leads to 



eo/2, where 

a = (1 -\- a/I -I- 4(72)/2. Going to the Cartan basis, the 
two generators which are common to these two SU{1, 1) 
algebras are 
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V, + eo 



and 
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where Tq is the generator of scale transformations. 



The choice of the third generator can be either Tj^ 



f = 



Ira Ei=i ^ i or formed from H as T\_ 
a E ».j=i T-^ — -IT • V,- where the f and i indices stand 
for "free" and "interacting" , respectively. Both of these 



satisfy the usual SU{1, 1) algebra: 

[r+,r_] = -2ro , [ro,T±] 



±T+ 



(9) 



The above generators obtained, after performing certain 
similarity transformations, are tailor made for the class 
of wave functions of Eq. ([6| in the scaling regime and 
are slightly different from [lO|. In particular, our algebra 
is valid in the infinite scattering limit, where the scaling 
law of Eq. ^ comes in to effect. The reason for the 
existence of two different ways of forming the algebra is 
the special property of the Euler operator Q = Ei ' 
in any dimension d, 

[Q,0'^] = rfO^ (10) 

which brings out the scaling dimension of a general opera- 
tor O. Thus, any two operators scaling with same degree 
can be combined to form a new algebra. For the case at 
hand, the two generators and T],, both have the same 
degree with respect to Tq. For instance, when the two 
body interaction is of the l/r^ type, scaling — > Afy, 
one can check that [To,v4] — showing the same 

behavior as 1/A^ Eili ^he fact that there exists a 
well defined Hamiltonian with good scaling properties as 
in ([?]), confirms that there are two different choices for 
5J7(1,1). 

A group theory view point to the zero energy bound 
states can be obtained by analyzing the ground state of 
the system, using the freedom of choosing a third gener- 
ator to perform a series of SU{1, 1) transformations. In 
this map, particular attention should be paid to the way 
in which the ground state and excited states transform 
under symmetry rotations. After an SU (1,1) rotation 
we obtain, 



He' 



H 



N 



(11) 



3 



where eo = + ^N{N — l)a. Putting the normahza- the free system obtained in this way obeys fractional 
tions back, remarkably, there is shift in the ground state 
energy, by ^fku. Thus, for an TV particle system, we see 
that the scaling exponent takes a universal formula inde- 



22|. A careful scrutiny of the energy reveals 



pendent of any details of the system as 7 = a ^^^^~^^ . 
The above analysis is in conformity with the sugges- 
tions in where 7 appears exactly in the actual shift 
of ground state energy of the N particle system, con- 
necting up well with the recent correspondence between 
non-relativistic conformal systems and their gravity du- 
als [1^ 17 1 . Each value of the scaling exponent is in 
gen eral related to the dimension of primary operators 
m,^- In this regard, the relation between the dimen- 
sion of the operators and the scaling exponent in our case 
is A = q; ^(^~i) ^ iK^ For two particles of species spin 
up and spin down, 7 can take values or —1 [5], which 
means a takes exactly the same values. For the case of 
three particles 7 « 0.22728 signifies the value is 0.07576 
for a. Since, both the SU{1, 1) algebras in Eq. Q can 
be embedded in the Schrodinger algebra, the exact map 
we constructed, confirms the duality between free and 
unitary fermions envisaged in [l3 | . The duality was orig- 
inally argued based on insights from the AdS/CFT cor- 
respondence [1^ regarding the existence of a pair of non 
relativistic conformal field theories with operators of dif- 
ferent dimensions. These features continue to hold for 
the Calogero type of models, as emphasized in . The 
duality between interacting and non-interacting systems 
is not restricted to the l/r^ potential, but to any gen- 
eral short-range potential which satisfies the scaling law 
in Eq. 

Now, owing to the reduction of [H,exp{—A/2}] to 
the commutator with respect to the Euler operator as 
Ei^i ■ Vi,exp{-4/2}] = ^exp{-l/2}, the operator 
T = $oexp{— A/2} diagonalizes the original Hamilto- 
nian H , giving 



statistics 

the possibility of exclusion statistics of the fermions near 
unitarity, which also points towards a mutual 1/r^ type 



short-range interaction 11 1 
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Finally, making another SU {1,1) rotation with respect 
to the generator T+ in Eq. one arrives at decoupled 
oscillators 



decoupled 



4e^? + Je-i + (^o-^^)- (13) 



This shows an exact map from the interacting model to 
the free model by a series of SU{1,1) transformations. 
Notice that only the ground state energy gets shifted and 
the excited states are same as the system of decoupled 
oscillators. In fact, in the absence of trap, the SU{1, 1) 
algebra can be used to make a unitary transformation, 
which takes particles interacting via the potential in ([3|) 
to free particles [2l|. Thus, U'^HU = 



where one takes U = e (^++^-)e^+^- . It is known that 



For example, in the three-body problem in a trap, in 
the unitarity limit, the interaction potential can actually 
be replaced by the following contact conditions 

V'(ri,r2,r3) = f — - -) A(R,,, r^) + 0(r,,) , (14) 

in the limit rij = jr^ — rj\ — > taken for fixed positions 
of the other particle k and of the center of mass Ry of 
i and j. In the molecule formation domain, A ^ and 
the wave function is singular, implying high probability. 
However, for the class of wave functions in Eq. which 
are of the Laughlin type, A = 0. Consequently, in the 
rij — > limit, the wave function has a dependence as in 
Eq. ([T]). The density of states in these cases is high, and 
the wave function actually vanishes when two particles 
come close to each other, much like the non-interacting 
system, signifying a dimerization, where the particles are 
weakly interacting. In the condensation regime the wave 
function should be singular, but yet normalizable, point- 
ing towards molecule formation. In this domain, the in- 
formation about the inter particle correlations is present 
in the Jastrow factor, and hence the statistics can come 
in to play. Universality of the degenerate gas on the other 
hand, suggests that the scaling exponent 7 does not de- 
pend on the specific statistics of the system. Neverthe- 
less, at unitarity, it might be useful to model the inter- 
action as coming from fractional exclusion statistics and 
intriguingly, the estimates for energy per particle are in 
good agreement with Monte Carlo simulations flll|. To- 
gether with the fact that these particles are interacting 
via a 1 /r^ potential, points towards a possible realization 
of exclusion statistics in this weakly interacting system of 
fermions ll|. As is widely known, exclusion statistics is 
realized in one-dimension [2^ by the Calogero model 24 1 
and in two dimensions, by the kinetic and potential en- 
ergy densities of fermions interacting with a zero-range 
potential in the mean- field approximation [25|. There is 
evidence emerging that such statistics can also be real- 
ized for cold atoms in three dimensions at unitarity [26| . 

The eigen functions of Hd are homogeneous and sym- 
metric functions, due to the presence of the Euler op- 
erator. The unitary irreducible representations of this 
SU{1, 1) algebra, with Tq diagonal can be constructed as 
well and point towards the existence of analytic scaling 
solutions, similar to the coherent states found in [27| . 
in a different context. The n-th energy level eigen 

function of H takes the form '^n — ^0 Pn , where, 
P„ = exp{— A/2}S'„ are totally symmetric, inhomoge- 
neous polynomials, with Sn denoting a choice of sym- 
metric polynomials. One particular choice can be made 



4 



by knowing that, and serve as the creation and an- 
nihilation operators respectively. The ground-state can 
be chosen: V^^o = , for i = 1, 2, • • • , A^. In two and 
higher dimensions, the excited states are built out of 
monomials depending on the square of the particle coor- 
dinates as Ylf {rf)^' , where ri; = 0, 1, 2, • • • . Here, is to 
be read as the hyperspherical coordinate which decouples 
from the rest of the system. The corresponding energy is 
E = 2J2ini+Eo leading to equally spaced excited states 
with universal frequency 2cj; the factor of 2 ensuing from 
the typical behavior. More over, these modes are 
independent of the initial N-particle system, with regard 
to the linear and non-linear response of atom-atom inter- 
actions in a trap ^ . The picture changes appreciably in 
one-dimension, where, the symmetric polynomials can be 
constructed without recourse to the hyperspherical coor- 
dinates. In particular, the excited states can be formed 
out of the particle coordinates themselves as Ylf {xi)"'' , 
which are understood to be symmetric with respect to 
the exchange of particle coordinates. Thus, the spec- 
trum of excited states turns out to be integrally spaced 
with frequency niuj^ for ni = 0, 1, 2, • • • . This leads 
us to predict that as opposed to higher dimensions, the 
resulting shift in the energies of the breathing modes in 
one-dimension are integral. This though mystifies the 
universality of the unitary fermi gas, suggesting a depen- 
dence of the breathing modes on the dimensionality of the 
trap, by singling out the special role of one-dimension. 
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